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The search for new interactions of neutrinos beyond those of the Standard Model may help to 
elucidate the mechanism responsible for neutrino masses. Here we combine existing accelerator 
neutrino data with restrictions coming from a recent atmospheric neutrino data analysis in order to 
f*"^ ' lift parameter degeneracies and improve limits on new interactions of muon neutrinos with quarks. 

In particular we re-consider the results of the NuTeV experiment in view of a new evaluation of 
its systematic uncertainties. We find that, although constraints for muon neutrinos are better than 



those applicable to tau or electron neutrinos, they lie at the few xlO 2 level, not as strong as 
previously believed. We briefly discuss prospects for further improvement. 

PACS numbers: 13.15.+g,12.90.+b,23.40.Bw 



(D . I. INTRODUCTION 

£■ 

The discovery of neutrino oscillations provides the only firm evidence for new physics we currently have [1], namely 
neutrino mass. This constitutes the most important discovery in particle physics in the last quarter century, implying 
. that the Standard Model, which correctly accounts for all other experimental data in particle physics, needs revision [2] . 
\ Precision is the keyword in neutrino physics today, including solar and atmospheric neutrino data analysis, as well as 
the determination of neutrino masses and mixing parameters at reactor and long-baseline accelerator experiments [3- 
5] . It is important not only to investigate the potential of current and upcoming long baseline experiments to determine 
the neutrino oscillation parameters, but also to probe for the possible existence of non-standard neutrino interactions, 
NSI, for short. The latter are expected in most models of neutrino mass generation, such as seesaw- type schemes [6] 
and will play a crucial role since they will shed light on the scale characterizing this so-far elusive mass generation 
, ^ . mechanism. The interest on probing for the existence of NSI has been growing over recent years, thanks also to the 
increasing precision of upcoming experiments [7, 8] and to the fact that the current oscillation interpretation is not 
yet fully robust [9] so one needs to scrutinize the interplay between oscillations and NSI in future experiments [10-13]. 
A wide class of non-standard neutrino interactions may be parametrized at low energies by the effective Lagrangian 

-£n f l = E s f a P p 2V2G F {D alp Lvp){fYPf) , (1) 

a/3 

where Gf is the Fermi constant and the parameters characterize the strength of the NSI. For simplicity, these 
are assumed to be real. The chiral projectors P denote {R, L = (1 ± 7 5 )/2}, while a and /3 label the three neutrino 
flavors, e, \i and r and / is a first generation charged SM fcrmion (e, u or d). 
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Collider experiments produce well-controlled and clean muon neutrino beams, with just a small component of 
electron neutrinos, while tau neutrino beams are unavailable. As a result one can expect that the muon neutrino NSI 
parameters, e p are better constrained in comparison with other neutrino flavors. This is indeed the case for the 
interaction with electrons [14, 15]. However, for the interaction of neutrinos with quarks, although neutrino nucleon 
scattering has a very long history since CDHS [16] and CHARM [17], the accuracy of these early experiments was 
worse than that of the more recent NuTeV experiment, which has measured the v^N interaction with a very high 
accuracy [18], reporting a discrepancy with the Standard Model predictions. While this may be interpreted as a hint 
of new physics, indicating at face- value a potentially non-zero value for the NSI parameters ej£ [19], uncertainties 
coming from QCD corrections might have been underestimated [20]. 

Moreover, as expected, limits on muon neutrino NSI coupling strengths [19] derived on the basis of 1-loop dressing of 
the neutrino effective four-fermion vertex can not be formulated rigorously. Indeed these are highly model dependent 
and, even when a full analysis is performed, including all the necessary diagrams required in order to obtain gauge- 
invariance, one does not strongly constrain the flavor changing parameters e q ^.' R and e q ^' R [21]. 

Recently, two new determinations of the Standard Model electroweak mixing parameter sin 2 9w from the NuTeV 
measurement have been presented [22, 23]. Given that the most stringent model-independent bounds on e^ a come 
from the NuTeV data, it is necessary to consider the effect of including these new corrections on the NuTeV results. 
Therefore, the main goal of this paper is to the review the current status of the constraints for the non-standard 
neutrino interactions in view of the larger uncertainties indicated by these recent papers. In so-doing we will not only 
take into account the new results mentioned above, but also combine the laboratory constraints with the restrictions 
inferred from the analysis of atmospheric neutrino data. In order to obtain these constraints we adopt as a simplifying 
working hypotesis that all NSI parameters other than e qR a vanish. We will see that restricting ourselves to a two- 
generation NSI global analysis we obtain relatively stringent constraints on the NSI interactions for muon neutrinos, 
at the few x 10~ 2 level, thanks largely to the interplay of atmospheric data. 

The paper is planned as follows. In the next section we will make a brief description of the neutrino nucleon 
scattering parameters that are relevant for our analysis. In section III, we will discuss the NuTeV data and their 
implications for the NSI parameters, while in sec. IV we discuss the previous experiments CHARM and CDHS. In 
section V we combine this information with constraints coming from the analysis of atmospheric neutrinos in order 
to obtain global constraints on the NSI parameters. Finally, in section VI we will give our conclusions and discuss 
the prospects for further improvement in the determination of NSI parameters. 



II. 9w MEASUREMENTS IN NEUTRINO— NUCLEON SCATTERING EXPERIMENTS 

Neutrino scattering experiments provide one of the most precise probes of the weak neutral current, and have been 
often been used to measure the electroweak mixing angle sin 2 9w- In particular, it has been shown that experiments 
with an isoscalar target are particularly convenient, since in this case the uncertainties due to unknown corrections to 
the QCD parton model cancel to a large extent [24]. For an isoscalar target of up and down-type quarks, the largest 
contributions to the neutral and charged current cross sections are related by isospin invariance and their ratio can 
be written as: 



where one introduces the ratio 



ajD^N -» M +A) 
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and the coupling constants arc denned as 

C<) 2 = (gf ? + (gf) 2 (5) 

with P=L,R. In the presence of NSI of muon neutrinos with quarks, these coupling constants are replaced by: 

(^) 2 = {gf + ^f + {gf + ^) 2 + El £ ^l 2 + El e ^l 2 ( fi ) 
{g«? = (gf + O" + Cflf + O 2 + E l^l 2 + E l<Sl a ( 7 ) 

The quantities R" and i2" have been measured in the past at the CDHS [16] and CHARM experiments [17]. 
Another well-known observable for the study of deep inelastic neutrino scattering on an isoscalar target is the Paschos- 
Wolfenstein (PW) ratio defined as [25] 

RPW = *{ Vtl N->n-X)-*{i>,N->^X) ~ 1-r V»> {9 »> ■ (8) 

This ratio is particularly useful because it depends very weakly on the hadronic structure of the nucleus target and it is 
largely insensitive to the uncertainties resulting from charm production as well as charm and strange sea distributions. 
However, the simultaneous measurement of neutral current cross sections for neutrinos and antineutrinos requires the 
use of separate neutrino and antineutrino beams. The NuTeV Collaboration makes use of this observable in order 
to measure the electroweak mixing angle. Actually, they measure experimentally the ratios R u and R p , shown in 
Eqs. (2) and (3), which later transform in the Paschos-Wolfenstein ratio Rpw through their fit procedure. 



III. CONSTRAINTS FROM NUTEV DATA 

The NuTeV collaboration used high statistics neutrino and antineutrino beams to measure their neutral and charged 
current cross sections on an iron target. Using a statistical separation of the NC and CC event candidates, based 
on the length of each event in the detector, NuTeV reported experimental values for R v and as mentioned in 
the previous section. A numerical fit of these results provides a measurement of the left and right-handed neutral 
couplings to the light quarks [18]: 

( g L) 2 = 0.30005 ± 0.00137 (9) 
{g^f = 0.03076 ±0.00110 (10) 

Notice the discrepancy with respect to the Standard Model expectations [1]: 

(#£)sm = 0.30399 ± 0.00017 {g*)l M = 0.03001 ± 0.00002 , (11) 

due to the NuTeV preference for a lower effective left-handed coupling, almost 3er away from the SM expectation. 
Similarly, a value for the electroweak mixing angle in the on-shell scheme was obtained from a numerical fit: 

sin 2 6 W = 0.22773 ± 0.00135 s * Qt ± 0.00093^* (12) 

which is 3cr away from the value determined in global precision electroweak fits [1]: 

sin 2 6 W = 0.22292 ± 0.00082 (13) 



Ascribing the discrepancy between the NuTeV result for the left coupling (g^) 2 and its prediction within the 
Standard Model to the existence of a non-zero NSI-opcrator, the authors in Rcf. [19] obtained a positive hint for 
non-zero values of the left flavor-conserving NSI couplings e d J^ and e^. On the other hand, for the case of flavor 
changing NSI parameters they obtained limits on the NSI |e^ ;fl i couplings. 
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FIG. 1: Current status of sin 2 9w determinations from Refs. [1, 18, 22, 23]. 



However, several corrections and theoretical uncertainties coming, for instance, from nuclear effects and next-to- 
leading-order corrections were neglected in the original NuTeV Collaboration analysis. Although the presence of NSI 
could explain the discrepancy of the NuTeV results with the SM, before claiming that they provide a hint of new 
physics, one needs make sure that all uncertainties are carefully taken into account. As a matter of fact several 
attempts have been made to interpret the NuTeV results just in terms of conventional physics, see for example 
Refs. [20, 26]. 

Here we re-analyze the impact of taking into account carefully the uncertainties in deriving restrictions on physics 
beyond the SM as parametrized by the NSI Lagrangian in Eq. (1). In particular, we will focus on the two most recent 
re-analysis of the NuTeV data given in Refs. [22, 23] in order to obtain constraints on neutrino NSI coupling strengths. 

The NNPDF Collaboration reports more precise estimates of the strange and anti-strange parton distribution 
functions, leading to a new value for the electroweak mixing angle [22]: 



sin 



= 0.2263 ± 0.0014 stat ± 0.0009 sys ± 0.0107 PjDFs 



(14) 



On the other hand the analysis in Rcf. [23] takes into account three different corrections coming from nuclear effects, 
due to the excess of neutrons in iron, charge symmetry violation arising from up and down quark mass differences, and 
strange quarks. With all these corrections included, the following value for the electroweak mixing angle is extracted: 



sin 2 6 W = 0.2232 ± 0.0013 stat ± 0.0024 s 



(15) 



One can see from Eqs. (14) and (15) that the two rc-calculated estimates for the electroweak mixing angle, with 
larger uncertainties, arc consistent with the SM prediction, as seen in Fig. 1. In order to determine the resulting 
restrictions set by NuTeV data on the strength of NSI couplings, in view of the corrections discussed above [22, 23], 
we adopt the Paschos-Wolfenstein ratio Rpw- We perform a \ 2 analysis using the Paschos-Wolfcnstcin ratio derived 
from each re- analysis of NuTeV data, and we compare these results with the Standard Model prediction for Rpw ■ 
As a result of this simple statistical analysis, and allowing for one non-zero NSI coupling at a time, we obtain the 
following constraints at 90% C.L. : 



-0.017 < & < 0.025 



0.84 < ei L u < 0.88, 



-0.72 < < -0.67 & -0.031 < < 0.020, 



-0.24 < e d * < 0.088. 
-0.058 < < 0.063 



(16) 



< < 0.36. (17) 



when using the results in NNPDF, and 



-0.005 < et L u < 0.005 
-0.71 < < 0.70 & 



0.86 < eft < 0.87, 
0.006 < < 0.006, 



-0.17 < e d * < 



-0.11 & 



-0.014 < 



0.042 < < 0.025.(18) 
uR < 0.016 & 0.28 < < 0.31. (19) 



/'/' 
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FIG. 2: Allowed regions at 90% C. L. and 3<r for the vectorial and axial NSI couplings of neutrinos with down (left panel) 
and up-type (right panel) quarks. Colored regions are obtained using NNPDF corrections [22], while empty lines correspond 
to Bentz et al [23]. 

for the Bentz et al. case. For both analyses we observe two allowed regions for most of the NSI couplings, reflecting 
the degeneracy seen in Fig. 2. 

Allowing now for the simultaneous presence of left and right-handed NSI neutrino couplings we obtain the results 
given in Fig. 2. For convenience, we present the constraints in terms of the vector and axial couplings instead of using 
the L, R basis. This makes it easier to proceed when combining with the atmospheric data (see below). The left panel 
holds for the case of a down-type quark, while the right one holds for an up-type quark. One can see that, in both 
cases, there is a two- fold degeneracy in the values of the NSI parameters, which is easily understood from Eq. (8) 1 . 
On the other hand, for the case of flavor changing NSI, our analysis gives the allowed region shown in Fig. 3. In this 
case the allowed regions for the two chiral components of the flavor-changing NSI couplings of u- and d-type quarks 
are the same, hence we show our results in only one plot. Note that here the allowed region is a hyperbola centered 
at the origin (ejX , £^t) = (0,0), displaying a remnant of the two-fold degeneracy seen for the flavor conserving case. 
In agreement with Eqs. (6)-(8), our analysis constrains the product e q ^ T ■ e^, and therefore, one coupling could be 
of order one, provided that the other one is small enough. As before, the different widths for the allowed regions for 
each re-analysis of NuTeV data reflect the experimental errors calculated in each case. 

IV. CONSTRAINTS FROM CHARM AND CDHS 

Both the CHARM and CDHS experiments have measured semileptonic neutrino scattering cross sections. Although 
these early experiments have a lower sensitivity than that reached at the NuTeV experiment, their data arc still useful 
in our global analysis. While these two experiments can not compete with the NuTeV sensitivity in constraining the 



1 Note that in terms of the NSI couplings eL, eL, this equation represents a hyperbola centered at (e^, = (~g% > -gfi )■ After the 

change of variables to (e*^, the hyperbola is rotated to the one shown in Fig. 2. 
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FIG. 3: Same as Fig. 2 for the flavor changing couplings of up- and down-type quarks with neutrinos. Colored regions are 
obtained using NNPDF corrections [22], while empty lines correspond to the Bentz et al. re-analysis [23]. The small region in 
the colored ellipse results from the analysis of CHARM and CDHS results, see discussion in Sees. IV and VI. 



Experiment 


Observable 


measurement 


SM prediction 


CDHS [16] 


R v 


0.3072 ± 0.0033 


0.3208 




ir 


0.382 ±0.016 


0.381 




V 


0.393 ± 0.014 




CHARM [17] 


RY 


0.3093 ± 0.0031 


0.3226 




R? 


0.390 ±0.014 


0.371 




r 


0.456 ±0.011 





TABLE I: Ratios of neutral to charged currents measured by CHARM and CDHS compared with the SM prediction. 



strength of non-universal NSI couplings, they play an important role in restricting that of the flavor changing NSI 2 
and for this reason we consider them in our analysis. The measurements of these experiments are summarized in 
Table I. Using the information given in this table we perform a \ 2 analysis, including the correlation between the 
neutrino and antineutrino ratios of neutral to charged currents given by the parameter r. Our \ 2 function in this case 
is given by 

X 2 = E^= ~ Knsi)(^(^ ~ Ri,NSi) (20) 

where i runs for either CHARM or CDHS and j, k stands for R v and R u . The results of our computation are shown 
in Fig. 3. Although the inclusion of these data strongly constrains the axial flavor changing NSI parameters, the 
restriction follows from a discrepancy between the experimental and the theoretical value of R v , as can be seen from 



2 The only caveat, however, is the bad quality of the resulting fit, given that these old experiments are not in good agreement with the 
SM, and that the addition of NSI only makes it worse. 



7 



Table I. As already mentioned, this results in a very a poor fit, hence we warn that constraints obtained for this 
case should be considered as less robust than those obtained for non-universal NSI discussed above. For the sake of 
completeness, however, and in analogy with Eqs.(16)-(19) we quote the constraints that would be obtained at 90% 
C.L. by combining CHARM/CDHS with the re-analysis of NuTcV data performed by the NNPDF Collaboration (or 
Bentz et al.): 

-0.023(-0.023) < ef T < 0.023(0.023), -0.039(-0.036) < ef T < 0.039(0.036). (21) 

V. COMBINING WITH ATMOSPHERIC DATA 

Atmospheric neutrino data are well described by the standard mechanism of neutrino oscillations [4, 5]. As a 
result, if neutrino NSI with matter exist, they must at best play a sub-leading role in the description of atmospheric 
neutrino data. We now discuss atmospheric neutrino propagation in the presence of non-standard neutrino-matter 
interactions, containing both flavor-changing and flavor-diagonal components [27-29]. In the context of such hybrid 
scheme, the presence of neutrino NSI would affect atmospheric neutrino propagation in matter introducing an extra 
term in the evolution Hamiltonian. In the simple two-neutrino description the NSI contribution to the Hamiltonian 
will be given by 

/ £ iv g v \ 

H NSI = V2G F N q mm £ (22) 

with q=u, d and = + e^, where one notes that neutrino propagation is only sensitive to the vectorial NSI 
couplings. This provides important information complementary to that coming from the accelerator experiments, see 
Eqs. (6) and (7). 

The Super-Kamiokande neutrino data have been analyzed under this assumption in Refs. [27-29]. Up to now no 
evidence of NSI has been found in the atmospheric data sample and, as a result, one gets upper bounds on the 
magnitude of the NSI coupling strengths. Here we will include in our analysis the most recent published results using 
the full atmospheric SK-I and SK-II data sample [29] which leads to the following bounds: 

0.007 < ef r < 0.007 (23) 
\e%-e%\ < 0.042, (24) 

at 90% C.L. (1 d.o.f.). These limits have been obtained for the case of neutrino NSI with down-type quarks. Taking 
into account the chemical composition of the Earth in the PREM model [30] , we can rewrite the bounds on Eq. (24) 
in terms of the neutrino NSI couplings with up- type quarks to a good approximation. According to Ref. [31], the ratio 
between the density of down and up quarks Nd/N u takes a value of 1.009 (1.047) in the mantle (core) of the Earth. 
Using the average value: Nd/N u = 1.028 one obtains the following 90% C.L. bounds on the neutrino - up quark NSI 
couplings: 

0.007 < ef T < 0.007 (25) 
\e%-e%\ < 0.043 (26) 

We now go one step further combining also these results with those of the previous accelerator analysis. The 
constraints derived from such a combined analysis are shown in Figs. 4 and 5 for the case of flavor-diagonal and 
flavor-changing NSI, respectively. As before, for the non-universal NSI couplings we consider NuTeV, while for the 
flavor-changing couplings we also include the results of CHARM and CDHS. As seen in the previous section, for the 
case of flavor changing NSI the NuTcV accelerator results and atmospheric neutrino bounds are nearly the same for 
up- or down-type quarks, and therefore we only present the plot for down-type quarks. One sees that the atmospheric 
neutrino data provide quite a useful tool to constrain the vector NSI coupling parameter, allowing to remove the 
degeneracy discussed before. 
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FIG. 4: Allowed regions at 90% C. L. and 3<r for flavor-diagonal NSI couplings of down- (left panel) and up-type quarks (right 
panel) with neutrinos obtained by combining accelerator data with results from the atmospheric data analysis. We use the 
same conventions as in Fig 2. 



Once we make the projection of the two parameter analysis into one parameter, we obtain the following constraints 
at 90% C.L. when considering the reanalysis of the NuTcV anomaly reported by Rcf. [22] ([23]): 

-0.042(-0.042) < e<w < 0.042(0.042), -0.091(-0.072) < eft < 0.091(0.057). (27) 

-0.044(-0.044) < e% < -0.044(0.044), -0.15(-0.094) < eft < 0.18(0.14). (28) 

For the case of the flavor changing parameters the results are nearly the same for the two different re-analysis 
of NuTeV data, since the determination of the vectorial couplings eft are dominated by the atmospheric neutrino 
analysis, while the axial couplings eft are mainly given by the accelerator experiments CHARM and CDHS. One 
finds, 

-0.007 < eft < 0.007, -0.039 < eft < 0.039. (29) 
-0.007 < eft < 0.007, -0.039 < eft < 0.039. (30) 

A more complete full-fledged three-neutrino global analysis in the context of neutrino oscillations plus non-standard 
interactions will introduce three additional NSI couplings, namely e ee , e e)Jl and e eT and, as a result, the bounds would 
be correspondingly weaker. 

VI. DISCUSSION AND FUTURE PROSPECTS 

We have re-analysed the constraints on novel non-standard neutrino interactions of muon neutrinos with quarks 
which follow from current accelerator experiments. In particular, we have re-considered the results of the NuTcV 
experiment in view of a new evaluation of the NuTeV systematic uncertainties. We have combined the restrictions 
following from accelerator data with those coming from recent atmospheric data analysis, which plays a crucial role in 
removing degeneracies. We have found that, although constraints for muon neutrinos are better than those applicable 
to tau or electron neutrinos, they are not as strong as previously believed. While previous authors report bounds of 
the order of few xlO -3 , we find that, taking into account the current estimates of the NuTeV uncertainties, constraints 
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0.02 



FIG. 5: Same as Fig. 4 for the flavor changing NSI couplings with down-type quarks, with nearly the same results for up-type 
quarks. The two re-analyses on NuTeV data considered here lead to the same results, since the analysis is dominated by 
atmospheric and CHARM + CDHS data. 



Global with NuTeV reanalysis 


NSI with down 


NSI with up 




NU 




NU 


NNPDF [22] 


-0.042 < e% 
-0.091 < eft 


< 0.042 

< 0.091 


-0.044 < e% < -0.044 
-0.15 < < 0.18 


Bentz at al. [23] 


-0.042 < eft 
-0.072 < eft 


< 0.042 

< 0.057 


-0.044 < e% < -0.044 
-0.094 < < 0.14 




FC 




FC 


NNPDF/Bentz et al. 


-0.007 < e% 
-0.039 < eft 


< 0.007 

< 0.039 


-0.007 < e$ < 0.007 
-0.039 < < 0.039 



TABLE II: Allowed values for NSI parameters from a global analysis of NuTeV data using NNPDF and Bentz et al. corrections, 
combined with atmospheric data. The top of the table corresponds to non-universal (NU) NSI parameters and the bottom to 
flavor-changing (FC) NSI. In the latter case CDHS and CHARM data are also included. 

are now of the order of few xlO -2 . Our results are summarized in Table II. Even weaker constraints would hold in a 
generalized three-generation framework. Although muon neutrino NSI couplings are often neglected [32], one should 
keep in mind the weakness of the existing limits. 

Note that these results will not be improved by the inclusion of data from the long baseline experiment MINOS [33], 
given its poor sensitivity to matter effects when compared with that of atmospheric data, at least in a two-neutrino 
analysis. As a result MINOS will not significantly restrict the neutral current NSI couplings discussed here, and the 
same is expected for OPERA [34]. However, future "clean" measurements of the oscillation parameters combined 
with the NSI-contaminated parameters determined by atmospheric experiments would be helpful. Certainly it would 
also help considerably to have a very long baseline setup as envisaged in the Long Baseline Neutrino Experiment 
(LBNE) project currently proposed, capable of probing Earth matter effects with enhanced sensitivity [35]. 

Last, but not least, we note that a new high statistics neutrino scattering experiment, NuSOnG has been pro- 
posed [36], with the same goals as NuTeV. One of its main tasks would be to probe new physics in the neutrino 
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couplings. Expectations are that NuSOnG would improve the NuTeV sensitivities on muon-neutrino quark scattering 
by a factor 2 at least, which would translate into a similar improvement on the NSI sensitivities with respect to those 
obtained from NuTcV. 
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